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Ahstract 


The  Melnikov  theory  of  perturbations  of  Hamiltonian  systems  con¬ 
taining  homoclinic  orbits  is  extended  to  systems  containing  canonical  vari¬ 
ables  belonging  to  the  coadjoint  orbits  of  a  Lie  group.  This  is  applied 
to  the  free  rigid  body  with  attachments  and  to  the  nearly  symmetric  top. 
These  syj‘ems  are  thereby  shown  to  have  transverse  homoclinic  manifolds 
in  an  appropriate  return  map  and  therefore  have  complex  dynamics.  In 
particular,  the  heavy  top  and  rigid  body  with  one  attachment  are  shown  to 
contain  horseshoes  and  therefore  have  no  additional  analytic  integrals, 
while  the  rigid  body  with  several  attachments  exhibits  Arnold  diffusion. 


1.1 


§1 .  Introduction 

This  paper  presents  theorems  which  establish  the  existence  of 
horseshoes  and  Arnold  diffusion  for  nearly  integrable  Hamiltonian  systems 
associated  with  Lie  groups.  The  methods  are  based  on  our  two  previous 
papers,  Holmes  and  Marsden  [1981b,c].  The  two  main  examples  treated  here 
are  as  follows: 

1.  The  rigid  body  with  attachments.  This  system  has  horseshoes 
(withone  attachment)  and  Arnold  diffusion  (with  two  or  more 
attachments) . 

2.  A  rigid  body  under  gravity,  close  to  a  symmetric  (Lagrange)  top. 

This  system  is  shown  to  have  horseshoes  (and  hence  is  not  integrable). 

The  main  new  feature  present  here  is  the  presence  of  Lie  groups.  Both 
the  symmetry  groups  and  the  basic  phase  spaces  involve  Lie  groups  and  our 
perturbation  methods  must  be  modified  to  take  this  into  account.  As  in  our 
previous  work,  the  results  hinge  on  reduction  together  with  a  method  of 
Melnikov.  This  is  used  to  analyze  the  perturbation  of  a  homoclinic  orbit  in 
an  integrable  Hamiltonian  system.  In  example  one,  the  unperturbed  system 
is  the  free  rigid  body  which  has  a  homoclinic  orbit  lying  on  a  sphere.  This 
sphere  arises  as  the  coadjoint  orbit  for  the  rotation  group  S0(3)  and  the 
computation  of  Poisson  brackets  needed  in  the  Melnikov  theory  is  most  easily 
done  using  the  KKS  (Kostar.t,  Kirillov-Souriau)  theory  of  coadjolnt  orbits. 

This  theory  is  reviewed  in  Section  2. 

After  removing  the  angular  momentum  about  the  vertical  axis,  the  cor- 

*  2 

rect  phase  space  for  a  rigid  body  under  gravity  is  T  S  the  cotangent  bundle 
of  a  sphere.  This  and  its  connection  with  Euler  angles  and  coadjolnt  orbits 
in  the  Euclidean  group  is  explained  in  Section  3.  This  section  thus  sets 


1.2 


up  the  basic  phase  spaces  needed  In  the  analysis  of  our  second  example. 

Section  4  develops  the  Melnikov  theory  when  the  phase  space  is  a 
product  of  a  KKS  phase  space  and  a  set  ov  action  angle  variables.  This  is 
applied  to  the  rigid  body  with  attachments  in  Section  5. 

Section  6  develops  the  Melnikov  theory  for  systems  on  a  phase  space 
where  the  unperturbed  system  admits  an  symmetry  and  has  a  homoclinic 
orbit  in  the  reduced  phase  space.  This  generalization  does  not  assume  the 
original  phase  space  is  simply  a  product.  This  generalization  is  needed 
for  and  is  applied  to  thenearly  symnetric  heavy  top  in  Section  7. 

The  two  examples  presented  here  were  selected  because  of  their  physi¬ 
cal  interest  and  because  they  provide  good  models  for  how  the  Melnikov 
theory  must  be  modified  for  systems  with  topologically  nontrivial  phase  spaces. 

The  presence  of  horseshoes  in  the  motion  of  a  nearly  symmetry  heavy  top 
implies,  amongst  other  things,  that  the  dynamics  is  complex  and  cannot  be 
captured  by  averaging  methods  (cf.  Akulenko,  Leshchenko  and  Chernousko 
[1979]},  that  the  dynamics  has  periodic  orbits  of  arbitrarily  high  period 
embedded  in  an  invariant  Cantor  set  and  that  the  system  admits  no  additional 
analytic  integrals.  The  latter  fact  is  consistent  with  known  classical  re¬ 
sults,  but  the  existence  of  horseshoes  is  a  stronger  and,  we  think,  more  sig¬ 
nificant  assertion.  Ziglin  [1981]  recently  proved  a  general  nonintegrability 
theorem  that  includes  the  non-standard  rigid  body  cases  but  his  proof  pro¬ 
ceeds  along  different  lines  and  does  not  seem  to  yield  horseshoes. 

We  expect  that  dissipative  and  forcing  terms  added  to  these  systems 
can  be  dealt  with  along  the  lines  of  Holmes  and  Marsden  [1981b]. 

Acknowledgements.  We  thank  Allan  Kaufman,  Steve  Smale  and  Alan  Weinstein 
for  motivation  and  several  useful  comments. 


2.1 


* 


Since  the  basic  paper  of  Arnold  [1966]^ Lie  groups  have  played  an  im¬ 
portant  role  in  the  construction  of  phase  spaces  and  the  symmetry  properties 
of  some  important  mechanical  systems.  For  systems  such  as  the  rigid  body 
one  wishes  to  realize  the  classical  Euler  equations  as  Hamiltonian  equations 
on  an  appropriate  phase  space  to  apply  Hamiltonian  perturbation  techniques. 

For  a  rigid  body  free  to  rotate  ahout  its  center  of  mass,  the  basic  phase 
space  is  T  SO (3)  which  is  conveniently  parametrized  by  the  Euler  angles 
(4>,<M)  and  the  corresponding  conjugate  momenta  (p, ,  p,.  pQ).  However  the 

9  tp  U 

SO (3)  symmetry  of  the  problem  enables  one  to  reduce  this  3  degree  of  freedom 
problem  to  a  one  degree  of  freedom  system  whose  phase  space  is  a  sphere  in 
body-angular  momentum  space. 

For  purposes  of  this  paper  we  are  interested  mostly  in  the  Poisson 
bracket  structure  since  it  is  this  which  directly  enters  the  Melnikov  theory. 

In  what  follows  we  describe  the  Poisson  structure  associated  with  a  Lie  group 
and  its  connection  with  coad joint  orbits.  (For  the  same  theory  with  emphasis 
on  the  symplectic  geometry,  see  Abraham  and  Marsden  [1978,  Chapter  4].)  Proofs 
may  be  readily  supplied  by  referring  to  the  preceeding  reference,  to  Arnold 
[1978]  and  to  Guillemin  and  Sternberg  [1980]. 


(a)  Generalities  on  KKS  structures 

Let  G  be  a  Lie  group  and  0^  its  Lie  algebra.  For  £,n  edj.,  U.a] 

* 

denotes  the  Lie  bracket  of  £  and  n-  Let  denote  the  dual  space  of  crp. 
★  * 

For  F  :(}  •»  IR,  and  the  variable  in  oq.  denoted  by  p,  define 

OF(y).V  =  (V,  f£>  , 


(2.D 


2.2 


where  <,  >  denotes  the  pairing  between  and  Oyt  and  DF(y):  qj*  -* 


IR 


is  the  usual  (Frechet)  derivative.  It  is  understood  that  evaluated 

at  the  point  u. 

The  KKS-bracket  of  two  functions  F,  G  :oj,*  ■*  IR  is  defined  by 


UF,G}Kw)  ■  -  f] 


(2.2) 


This  bracket  makes  the  smooth  functions  from  CJ.*  to  IR  into  a  Lie  algebra. 
The  only  non-obvious  condition  is  Jacobi's  identity. 

Next  we  describe  the  relationship  between  the  KKS  bracket  and  coadjoint 
orbits.  For  g  €  G,  let  Adg  be  the  adjoint  representation  (the 

linearization  of  the  map  1^:  h  -*  ghg"^  at  h  =  identity)  and  Ad*_^  :<£.*  ■*■<>}+ 
the  coadjoint  representation.  For  y^  let  <2  -  (Ad*  -jUq  1  g  e  G  }  be 

the  orbit  of  yQ.  A  theorem  of  Kirillov,  itostant  and  Souriau  states  that  0 
is  a  symplectic  manifold.  We  now  describe  the  symplectic  structure.  For 
y  GO,  tangent  vectors  to  0  and  y  have  the  following  form:  let  ^ 
and  define  £(y)  e  qp  by  n  -"<0,  U»n]  >.  Then  £(y)  is  tangent  to  0  at 
u.  The  formula 


w^U(y),  n(y))  3  y»  U>n3  > 


(2.3) 


defines  a  symplectic  form  on  <?.  Formulas  (2.2)  and  (2.3)  are  related  as 


follows:  for  F,G  :  *  IR  , 


{(F,G) )  jO  -  (FlO  .  G| 0}  , 


(2.4) 


wh?re  {  is  the  Poisson  bracket  computed  from  the  symplectic  structure 

( 2 .  .1 )  on  0  . 

If  H  ■*  IR  is  a  given  Hamiltonian,  there  is  a  unique  vector  field 
*  * 

on  such  that  any  function  F  :  oy  -*■  IR  which  evolves  along  the  flow 


of  satisfies 


F  =  {{F,H}}  . 


(2.5) 


Explicit  equations  of  motion  for  y  can  be  computed  directly  from  (2.5) 

and  (2.2)  in  examples,  without  the  computation  of  coadjoint  orbits,  by  letting 
F  be  coordinate  functions  on  y*.  The  equations  so  obtained  are  thus 


y  3  XH(y) 


(2.6) 


The  vector  field  Xu  is  necessarily  tangent  to  the  coadjoint  orbits;  thus 
ri 

if  y(0)  e  0,  then  y(t)  €  0  as  well .  Furthermore,  X^jo  =  X^  where 
X^  is  the  Hamiltonian  vector  field  on  0  computed  from  H  using  the 
symplectic  structure  (2.3). 

As  we  shall  see  in  the  next  sections,  the  classical  Euler  equations 
for  a  rigid  body  (with  or  without  gravity)  can  be  expressed  in  the  form  (2.5). 


Remark.  In  Arnold  [1966]  and  Ebin  and  Marsden  [1970]  it  is  shown  that  the 

equations  of  an  incompressible  fluid  also  fall  into  this  class  using  the 

qrouo  V  of  volume  preserving  diffeomorphisms  of  space.  For  compressible 
vo  1 

flow,  the  appropriate  group  is  the  semi-direct  product  of  diffeomorphisms  and 
functions.  (The  reason  semi -direct  products  are  relevant,  is  explained  ab¬ 
stractly  below).  When  appropriately  coupled  to  the  electromagnetic  field, 
equations  (2.5)  also  include  the  equations  of  plasma  physics  and  nsagnetohydro- 
dynamics  (see  Marsden  and  Weinstein  [1981]). 

For  the  rigid  body  free  to  rotate  about  a  fixed  point,  the  basic  phase 
space  one  starts  with  is  T*S0(3),  irrespective  of  whether  or  not  gravity  is 
present.  In  the  absence  of  gravity,  reduction  by  S0( 3)  leads  naturally  to 


the  KKS  Poisson  structure  on  the  Lie  algebra  of  SQ(3).  However  when  gravity 
is  present,  reduction  by  leads  to  the  KKS  Poisson  structure  on  the  Lie 
algebra  of  the  Euclidean  group;  i.e.  the  semi-direct  product  S0(3)  €  IR  . 
Proofs  of  assertions  made  in  the  following  may  be  found  in  one  or  more  of 
Abraham  and  Marsden  [1978],  Guillenin  and  Sternberg  [1580]  and  Ratiu  and 
van  Mcerbeke  [1981]. 


(b)  KKS  Structures  for  Reduction  by  a  Subgroup’ 

■k  + 

Let  6  be  a  Lie  group  and  T  G  its  cotangent  bundle.  Let  v  €  cfr 
* 

and  let  G,,  =  [g  £  G|Ad  =  v}  be  the  isotropy  subgroup  of  v.  Now 

v  q"  ■ 

3  * 

G  acts  on  G  by  left  translation  and  hence  on  T  G.  This  action  has  an 


Ad  -equi variant  momentum  map 


J  :  T  G  '*•0^ 


where  is  the  Lie  algebra  of  G^.  In  fact,  for  a  €  T  G  and  i  €  0^ 


uu(cg).c>  -  <v  yg.«> 


(2.7) 


where  :G  ■*  6  denotes  right  translation  by  g  and  (  ,  >  denotes  the 

♦  i k 

pairings  between  <^v  and  <3^  or  T  G  and  T  6  as  appropriate. 

*  _  * 

Now  let  u  G  ^  and  u  be  its  restriction  to  We  shall 

henceforth  make  the  (generic!  assumption  that  G  is  abelian.  Then  the 

V  - - - 

* 

reduced  space  of  T  G  relative  to  the  action  of  G  is 


The  heavy  rigid  ton  is  done  two  separate  ways  in  Section  7,  namely  in 
terms  of  Euler  angles  and  in  terns  of  the  KKS  structure.  Those  who  wish 
only  to  read  the  Euler  angle  proof  may  omit  the  rest  of  this  section  and 
the  third  part  (c }  of  the  next  section. 


2.5 


PjtV*  (oVGiI/Sv  (2.8) 

* 

Let  <<,  >>  be  a  left  invariant  metric  on  G  and  K  :T  G  -*•  IR  the  correspond 

ing  kinetic  energy  function:  K(ag)  =  l(<agJag>>.  Let  be  the  one  form 

on  G  defined  at  g  e  G  by  minimizing  K  over  the  affine  space  of 

such  that  Jv(a  )  =  y. 

9 

Remark.  If  G^  =  S1  and  Z  the  element  corresponding  to  y  via 

<<,  >>  i.e.  *  uCn)  for  all  n  e^-v»  then  for  v€  T^G, 

yS)-v  ■  «TeRg5.  v». 

The  function  K(a  (g))  =  V  (g)  is  called  the  amended  potential, 
y  y 

The  map  a  -*•  «  -  a,  induces  a  symplectic  diffeomorphism  of  P- 
q  q  P  r  y»v 

with  T  (G/Gv)  by  Theorem  4.33  of  Abraham  and  Marsden  [1978]. 

* 

Remark.  In  general,  the  symplectic  form  on  T  (G/Gv)  is  the  canonical  one 
plus  the  "magnetic"  field  dc^.  For  the  examples  in  this  paper  du^  =  0. 
However  a  f  0  and  it  is  necessary  to  use  it  to  form  the  amended  potential. 
If  this  is  added  to  the  standard  potential  for  the  heavy  top,  we  recover 
the  effective  potential.  As  we  shall  see  below,  this  agrees  with  that  in 
standard  texts  (Goldstein  [1980],  formula  5-60,  p.  215)  and  is  a  special 
case  of  Theorem  4.5.6  of  Abraham  and  Marsden  [1978], 

Now  consider  the  semi -direct  product  G  of  G  with  the  additive 
with  G  acting  on  01  by  the  adjoint  action.  For  v  € 


op  ,  let 


2.6 


fv  :  G  x  oi  +  IR 


fV(g,C)  *  <V,  Adg_,c  >  (2.10) 

i.e.  f  (g* * )  =  Ad^v-  Write  f^(g)  =  fv(g,£).  From  the  identity 

fvLi,]9>t)  ■  fv(g,  Ad  5)  (2.11) 

1  91 
* 

we  see  that  G  x  acts  on  T  G  by  the  right  action 

{g,£Wh  -  t  v.(TRg)*.ah  (2.12) 


where  t  v  is  fiber  translation  in  T*G  by  the  differential  of  f^.  The 

,df£  *  5 

action  (2.12)  is  symplectic  and  has  an  Ad  -equi variant  momentum  map 


★  ★ 


A  :TS  + 


9  x 


given  by 


*  .  * 


A  <ag)  =  CCTeLg_1 )  ag,  Adg_.,v) 


(2.13) 


yc  fc 

We  shall  write  m  =  (T  l  , )  a  and  v  =  Ad  ,v.  (See  Abraham  and  Marsden 

e  g-i  g  g-i 

[1978],  Ex.  4.2c  and  Guillemin  and  Sternberg  [1980]).  We  note  that  the  ad¬ 
joint  action  of  G  x  on  $x  ?is 


Ad(g^)(n,c)  =  (Ad  ri,  Adg;  +  [C,  Adgn]) 


(2.14) 


and  that  the  Lie  bracket  on 


V  x 


[(^,n1)»  (C2,  n2)]  =  (CSp^l,  Urn2]  +  [nry)  (2.15) 


One  can  check  directly  from  (2.14)  and  (2.12)  that  Av  is  equivariant. 


(The  identity 


2.7 


AdgU,n]  =  [Adgc,n3 


is  useful  in  verifying  this). 


Observe  that  Av  is  invariant  under  the  left  action  of  6^;  i.e. 


for  h  e  G 


AV(TL*  ,  a  )  *  Av(a  ) 


(2.16) 


Thus,  A  induces  a  map 


One  readily  sees  that  F  is  a  diffeomorphism  of  P-  v  onto  the  orbit  ^  v* 


2.1  Theorem  A  :  P-  ->0  is  a  symplectic  diffeomorphism. 
y,v  y  ,v 


Remarks.  1.  This  result  is  due  to  Ratiu  [1981].  The  proof  we  give  is 
simpler,  being  motivated  by  Guillemin  and  Sternberg  [1980]. 

2.  For  the  Lagrange  top,  2.1  can  be  proved  by  a  direct,  but  messy, 
calculation  which  we  outline  in  the  next  section. 


Proof  of  Theorem  2.1.  Since  the  symplectic  form  on  P-  is  induced  from 

-  J  r  y,v 

* 

the  canonical  symplectic  structure  on  T  G  and  that  on  0  is  determined 

y  ,v 

by  the  KKS  structure,  it  suffices  to  show  that  Av  commutes  with  Poisson 
brackets.  This,  however,  is  a  general  fact  about  Ad*-equi variant  momentum 
maps  and  collective  Hamiltonians  proved  in  the  next  two  lemmas. 


'if  ^ 

■  be  an  Ad  equi variant  momentum  map  for, the  n- qht 


2.2  Lemma.  Let  J  :  P 


2.8 


Ihen  the  Hamiltonian  vector  field  Xpo j  for  the  (collective)  Hamiltonian 
F°  J  *.  P  -*•  IR  is  given  at  x  €  P  by 


where  e^is  evaluated  at  y  =  J(x)  and_  st< 
field  on  P  generated  by  the  Lie  algebra  element 


(2.17) 


Proof.  By  the  chain  rule  and  definition  of  the  functional  derivative, 


d(F°J)»v  =  dF(y)*dJ(x)*v 

A  A 


-'dJW-V 

where  vx  e  T^P  and  y  =  J(x)  as  above.  By  definition  of  a  momentum  map, 


<dj{x;.vx>|>=o,x((|)  v; 


where  w  is  the  symplectic  form.  Thus 


d(r  -  “4(l)p- 


which  means  X, 


[dy]. 


2.3  Lemma .  Under  the  assumptions  oc  the  proceeding  lemma,  if  F,G:^/  +  IR, 


{F°J.  G01}  =  {{F,G}}° J 


(2.18) 


Proof.  {F«0,  G°J}  =  co  fffi  ,  [£)  1 
-  xU6dJp  W;pj 


(2.19) 
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by  definition  of  Poisson  brackets  on  P  and  the  preceeding  lemma.  On  the 
other  hand,  the  KKS  bracket  is 

UF,G»<y)  =  c|£,  )  (2.20) 

On  the  other  hand,  Ad  equivariance  gives  (Abraham  and  Marsden  [1978,  Corol¬ 
lary  4.2.9]) 

<  J{x).  [§.  -  -t<J(x).|£>.  <J(x).  |£>>  (2.21) 

in  which  p  is  fixed  and  x  is  variable.  There  is  a  minus  sign  here  since 
d  is  generated  by  a  right  action.  Holding  fixed, 

d<JW-f)-vx  =  vx=<d'J(x)-vx.  f  > 

so  <d(x),  generates  the  same  Hamiltonian  vector  field  as  F°d,  as  in 
2.2.  Thus,  substitution  of  (2.21)  in  (2.20)  yields  (2.19).  ■ 

Remarks .  As  we  shall  see,  the  heavy  rigid  top  Hamiltonian  is  a  collective 

Hamiltonian  for  Jv.  This  is  compatible  with  the  reduction  picture:  a 

Hamiltonian  system  on  the  reduced  space  P-  can  be  written  in  the  form 

4  *  v 

F  = (2.22) 

where  H  is  written  in  terms  of  the  KKS  variables  (m,v)  . 


3.1 


§3.  The  rigid  body 

This  section  is  divided  into  three  parts.  The  first  part  explains 
how  to  write  the  equations  of  a  free  rigid  body  free  to  rotate  about  its 
center  of  mass  in  KKS  form  (2.5).  The  second  part  recalls  the  Euler  angle 
formulation  of  the  heavy  top  and  the  third  part  puts  it  into  KKS  form  giving 

4* 

the  explicit  relationships  with  Euler  angles. 


(a)  The  free  rigid  body 

* 

The  free  rigid  body  is  a  left  invariant  Hamiltonian  system  on  T  $0(3). 
where  S0(3)  is  the  group  of  proper  orthogonal  linear  transformations  of 

3 

IR  to  itself.  By  general  facts  about  reduction  we  know  that  the  equations 
of  motion  must  be  in  the  form  (2.5):  this  is  true  of  any  left  invariant 
Hamiltonian  system  on  a  Lie  group  G.  For  the  free  rigid  body  we  can  bypass 
the  Euler  angle  description  (the  relevant  formulas  are  given  in  Table  1 
below). 

The  Lie  algebra  so (3)  of  so (3)  consists  of  the  set  of  3  x  3  skew 

3 

symmetric  matrices.  We  identify  so (3)  with  IR  by  identifying 


v  =  (p,q,r)  e  1ST  with  v 


0  -r  q 

r  0  -p 

•q  p  0  J 


€  so(3)  (3.1) 


The  Lie  bracket  corresponds  to  the  cross  product  in  the  sense  that 

/v  /v  * 

[v,w]  =  (v  X  w)  . 


(3.2) 


The  third  part  may  be  omitted  if  desired  since  we  give  two  proofs  of  our 
main  result  for  the  heavy  top,  one  using  only  Euler  angles  and  the  other 
using  the  KKS  structure 


3.2 


where  {{  ,  }}  is  the  KKS  bracket.  In  the  present  case  this  bracket  becomes 

{{F,G}}(m}  =  -m«(7F  x  vG)  (3.7) 

for  F,G  functions  of  m. 

The  fact  that  the  equations  (3.6)  must  preserve  coadjoint  ov'fa 1 ts  amounts 
in  this  case  to  the  fact  that 


3.3 


lZ  =  m*  +  m|  +  m \  (3.8) 

is  an  (obvious)  constant  of  the  motion  for  (3.4).  In  terms  of  coadjoint 
orbits,  equations  (3.4)  are  Hamiltonian  on  each  sphere  in  m-space  with 

4* 

Hamiltonian  function  (3.5). 

The  flow  lines  are  given  by  intersectii.y  the  ellipsoids  H  =  constant 
with  the  spheres.  For  distinct  moments  of  inertia  the  flow  on  the  sphere 
has  saddle  points  at  (0,  ±A,  0)  and  centers  at  (±2.,  0,0),  (0,0, ±2,).  The 
saddles  are  connected  by  four  heteroclinic  orbits,  as  indicated  in  Figure  1. 


Figure  1.  The  spherical  phase  space  of  the  rigid  body  for  fixed  total 
angular  momentum  l  =  +  m^  +  m^;  1^  > 

TThe  coadjoint  orbits  are  spheres  because  the  coadjoint  action  of  SO (3)  on 
m  space  is  just  by  rotations.  The  KKS  symplectic  structure  is  proportional 
to  the  area  element  by  a  simple  calculation. 


3.4 


The  orbits  are,  of  course,  explicitly  known  in  terns  of  elliptic  func¬ 
tions.  The  orbits  of  the  most  Interest  to  us  are  the  heteroclinic  orbits  which 
are  given  as  follows.  The  four  orbits  lie  in  the  invariant  planes 


where 


m3  =  *  Jif  ml 


U  "  !3  *3  -  h 

a,  =  -i - ±  >o,  a0  =  <  0  and 

1  i2i3  2  W 


i  .  i 

a3 =  ~n —  >  °* 


3.1  Theorem.  The  heteroclinic  orbits  for  the  free  rigid  body  are  given  by 

+  /a 7  _ 

m^Ct)  =±ij~  sech(-/a^a^  Jit)  , 

mj(t)  -tl  tanh(-»^1a3  Jit)  ,  (3.9) 

m3(t)  =  +ji/~  sech(-/a1a3  Jit)  , 

/a7 

for  m3  *  +  m-j  and  by 


mj(t)  =  m|(-t),  m^t)  =  m^-t),  m"(t)  a 


-m3(-t) 


(3.10) 


for  m3  - 


f3 

'a1  V 


This  may  be  checked  by  direct  computation  or  by  consulting  one  of  the 
classical  texts. 


«*'*.**  si  ■Hra-'-wt 


3.5 


(b)  The  Heavy  Top:  Euler  Angle  Description 

We  now  recall  the  traditional  Euler  angle  description.of  the  heavy 
top  and  shall  locate  homoclinic  orbits  for  the  synwetric  (Lagrange)  top. 

The  Euler  angle  description  is  more  familiar  but  in  some  respects  the  KKS 
description  is  simpler.  For  this  reason  we  shall  present  both. 

Given  a  rotation  A  e  S0(3)  we  let  the  corresponding  Euler  angles  be 
denoted  (<j>,ij/,8)  using  the  conventions  of  Goldstein  [1980];  see  Figure 
2(a).  The  corresponding  conjugate  momenta  are  denoted  p^,  p^5  pQ  so  that 
U,<M»  P^.  P^,  P0)  coordinate  T  S0(3). 

We  let  m  denote  the  angular  momentum  in  the  body  and  let  v  =  A~^k 
where  k  is  the  unit  vector  along  the  spatial  z-axis.  We  assume  the  center 
of  mass  is  at  (0,0,5,)  when  A  is  the  identity.  The  vectors  (m,v)  are 
expressed  in  the  body  coordinate  system;  see  Figure  2(b). 

The  Hamiltonian  is 


3  m2 


H(m,v)  =  i  [  7*2-  +  Mgi.  v- 
2  j=l  3 


(3.11) 


where  M  is  the  total  mass.  When  written  in  terns  of  Euler  angles,  this 
becomes 

,2 


H 


1  ~  p<i'  cos  sin  ^  *  pe  sin  6  cos  jj* 


I 


p 

1^  sin  6 


((pe  -  cos  e)  cos  -  pe  sin  8  sin  \|>12  ?2 


I2  sin2e 


+  -r-  /  +  Mg£  cos  e 
*3  1 


(3.12) 


In  Table  1  below  we  summarize  the  relationships  between  ro,  vf  $, 

p,,  P,,  Pfl»  ^  and  0  for  convenient  reference,  but  in  this  sub¬ 
s'  v  a 

section  we  shall  take  (3.12),  or  equivalently  the  Lagrangian 


3.7 


l  =  -y  (<j>  sin  6  sin  i|>  +  0  cos  >j>)^  +  -j-  U  sin  0  cos  t|»  -  0  sin  4>)^ 

•  (  o 

+  -jjr  (<j>  cos  6+  tj>)  -  Mg*,  cos  0 


(3.13) 


as  our  starting  point,  and  work  exclusively  in  terns  of  Euler  angles. 

The  Hamiltonian  (3.12)  is  invariant  under  rotations  about  the  z-axis; 

i.e.  $  is  a  cyclic  variable,  so  is  a  constant  of  the  motion.  In 

other  terms,  the  momentum  map  for  this  action  is  p^,  p^,  pg)  = 

p^.  The  reduced  space  ^(p^/S1  is  parametrized  by  (t|>,0,  p  ,  pQ).  In 

★  2 

fact  this  reduced  space  is  identifiable  with  T  S  ,  the  cotangent  bundle  of 

the  two  sphere  on  which  (ip ,0 ,  p,,  p.)  are  canonical  coordinates.  In  fact, 

cs  we  shall  see  below,  Oj>,e)  are  a  system  of  spherical  coordinates  for  this 

sphere.  The  equations  of  motion  for  ip,  6  are  thus  just  Hamilton's  equations 

for  (3.12)  with  px  held  constant. 

<P 


Remarks  1.  The  reduction  here  is  in  accord  with  the  general  fact  that 
reducing  T*Q  by  S1  gives  T*(Q/S] ) ;  here  Q  »  $0(3)  and  S0(3)/S]  =  S2. 

2.  The  two  sphere  obtained  here  is  not  to  be  confused  with  the 
sphere  for  the  free  rigid  body  shown  in  Figure  1. 


For  the  symmetric  top,  1^  =  Ig,  <p  is  also  a  cyclic  variable  and 

p  is  constant  as  well.  In  fact  these  two  symmetries  commute,  so  we 
2 

have  a  torus,  T  ,  symmetry,  which  makes  the  system  completely  integrable. 

The  reduced  system  has  one  degree  of  freedom  in  (8,  p_),  The  reduced  space 

C7 

is  no  longer  a  manifold,  but  has  a  boundary  at  6  e  0,  n;  i.e.  it  is  T*[0,tj]. 


3.8 


These  singularities  In  the  reduced  space  correspond  to  the  fact  that  the 

4* 

level  set  p  a  constant  is  singular  at  0  *  0,ir  . 

With  =  I2,  (3.12)  becomes 


H  = 


-•5-  Up,  -  P,„  cos  6)2  +  p2  sin2e)  +  /  +  Mgjt  cos  e 


21^  s'm^e 


kvr4>  r'i>  —  re  i3 


2Il 


Mgs,  cos  6 


, (p* - r 9)2  A\ 

21-j  sin2e  *3 


(3.14) 


which  shows  the  amended  potential  explicitly. 


Remark.  In  these  coordinates,  the  abstract  formula  (2.9)  for  the  4, ^reduc¬ 
tions  gives  the  closed  one  forms  p^dp^  an<^  P^P^  respectively* 


3.2  Theorem.  Jf_  0  <  p  <  2vMgiI1  and  if  p  «  3  I^b,  the  reduced 

Lagrange  top  system  in  (9,  p0)  space  has  a  hyperbolic  saddle  points  at 

6  «  0,  p  =0  and  a  homoclinic  orbit  connecting  it  to  itself  given  by 
0  ■ 


,_2 


cos  9  =  1  -  y  sec’n 
L2 


IMS) 


where  0  0  2  Mgi/Ij  and  y  3  2-  -j. 


(3.15) 


Remarks  1.  A  top  with  0  <  p^  <  2,'MgS.^  for  which  the  vertically  spinning 
state  is  unstable  (a  saddle)  is  called  a  slow  top. 


fThis  is  a  special  case  of  a  general  fact  about  singularities  in  level  sets 
of  momentum  maps;  see  Arms,  Karsden  and  Honcrief  [1981]. 


)  stea 
atlon 


3.10 


2.  For  =  p^  note  that  the  potential  in  (3.14)  does  not  have  a 
singularity  at.  a  ■-  0,  nr,  so  these  end  points  do  not  cause  difficulties. 

3.  Note  that  0  <  y  <  2. 

The  theorem  is  most  easily  checked  by  using  the  energy  equation  (3.14); 
the  homoclinic  orbit  has  energy  H  =  MgJt  +  P^I^;  cf.  Goldstein  [1980], 
pp.  215,  216  and  Figure  3. 


n : 

I*  a 


(c)  The  heavy  top:  KK.S  description 

The  abstract  theory  in  the  preceeding  section  guarantees  that  the 
*  2 

reduced  space  T  s  for  the  heavy  top  is  symplectically  diffeomorphic  to  a 
coadjoint  orbit  in  the  semi-direct  product  SQ(3'  x  IR  ;  i.e.  in  the 
Euclidean  group  E.,.  The  Lie  algebra  is  denoted  e^. 

The  mapping  that  gives  this  diffeomorphism  is  just  the  map 


A  :  ($,  'M.  p^.  P^,  Pe)^  (m,v) 


(3.16) 


where  m  is  the  angular  momentum  in  the  body  and  v  is  the  orientation  of 
gravity  as  viewed  from  the  body.  Table  1  below  summarizes  the  explicit 
formulas  relating  these  quantities.  Table  2  summarizes  the  relationships 
between  the  "Euler  angle"  spaces  and  the  coadjoint  spaces. 

The  KKS  bracket  for  functions  of  (m,v)  is  given  by  (2.2)  and  (2.15), 
which  in  this  cas^  becomes 


{{F ,G}}(m,v)  =  -rn-(VmF  *  VmG)  -  v.(vmF  x  v„G  +  VvF  x  VmG) 


m 


m 


m 


(3.17) 


The  assertions  of  the  general  theory  can  be  checked  by  hand  in  this 
case.  F<rst  of  all,  it  is  easily  seen  that  the  KKS  bracket  equation 
F={{F.,H}}  with  H  given  by  (3.11)  and  the  bracket  by  (3.17)  yield  the 
equations 


W  1 


Table  1 


Formulas  Relating  Euler  Variables  and  KKS  Variables  for  the  Heavy  Top 


m-|  =  [(p  -p^  cos  e)  sin  ip  +  pQ  sin  e  cos  t^j/sin  0  3  I-jCI  sin  e  sin  ip  +  e  cos  ip) 

m2  =  COyP^  cos  0)  cos  ip  -  Pq  sin  0  sin  t|>]/sin  0  3  sin  ©  cos  ip  ~  6  sin  tp) 

♦  • 
m3  3  3  I3(.<J>  cos  6  +  \p) 

v.j  =  sin  e  sin  ip 

v2  =  sin  0  cos  ip 

v3  3  cos  0 


•  •  •  « 

=  m«v  3  I.j(<j>  sin  0  sin  \p  +  6  cos  tp)  sin  0  sin  \p+I^(<P  sin  0cos  tp  -0  sirup)  sin  ecos  ip 

*  m3  3  I3(<j>  cos  0  +  ip) 

pQ  =  (v2mi'’vim2^y^’v2  =  *1  ^  sin  9  s*n  ^  +9  cos  ^)cos  V  '  sin9  C0S1^  "9  sin  ip)  sin  \p 


•  i  m.  v  -  -l  tn^  v  t, 

x  _  t  I  I  j.  i  d  c 
<f>  -  f - +  f - 

in  7 - n  if 


m3  m3ml V1  m3m2V2 


‘3  i, 


e  » 


mlv2 


Vl 


M  /V1  2  l2  /V1  2 


3.12 


3.13 


a  a^rr^  -  Mg£  v2 
m2  =  a2mlm3  +  V1 


*  a3m}m2 


_  m3vi  m2v3 

1  ’  17 "  V 


v2  * 


ml v3  m3vl 


*1  *3 

m2',l  m1  v2 


I 


1 


(3.18) 


(3.19) 


A  slighlty  tedious,  though  straightforward  computation  shows  that  these  are 
equivalent  to  the  Hamiltonian  equations  for  (3.12).  In  fact,  if  F,  G  are 
functions  of  (m,v)  and  (m,v)  are  then  written  out  in  terms  of  the  Euler 
variables  then  an  explicit  (but  tedious)  computation  shows  that 


{{F,G}}(m,v)  =  {F,G>0^  (3.20) 

where  (  ,  ,  denotes  the  bracket  computed  for  the  e,i//  variables  (i.e. 

0  ,1 {I 

holding  p^  constant).  This  is  the  content  of  Theorem  2.1  for  this  special 
case. 

The  equations  (3.18)  and  (3.19)  have  Ivl  and  m.v  as  constants  of  the 

motion.  This  just  reflects  (a)  the  conservation  law  pA  *  constant  and 

<J> 

(b)  the  preservation  of  the  coadjolnt  orbit  by  the  KKS  equations.  The  condi¬ 
tions  llvll  =  1  and  m«v  ■  p  »  constant  also  give' the  explicit  Identification 

y 

*  ? 

of  the  coadjolnt  orbit  with  T  S  .  Indeed,  Ivl  ■  1  describes  the  unit 

sphere  S  and  m*v  ■  p  specifies  m  as  a  linear  functional  on  the  unit 
2  2 

normal  to  S  ,  leaving  m  restricted  to  TyS  free.  Thus  m  determines, 


■>:<  rsTvTJjfnxv'i 


v, T J « W: ^ rrP&*f7Fr - w y/y^v.  srv  r-.’tsww^' 


3,14 


*  2 

by  restriction,  an  element  of  TvS  .  Finally  observe  that  the  equations 

2 

for  v  in  Table  1  show  that  e»ip  are  spherical -type  coordinates  on  S 
0j>  has  been  rotated  by  rr/2  from  standard  conventions  on  spherical  co¬ 
ordinates). 

Finally  we  discuss  the  Lagrange  top  in  the  KKS  picture.  For  I-j  =  I2 
the  invariance  is  rotations  about  the  3  axis.  Tnis  action  corresponds 
to  the  action  of  rotation  through  tp  in  the  Euler  angle  picture,  as  is 
easily  seen.  Also,  the  momentum  map  can  be  directly  checked  in  the  KKS 
picture  to  he  just  m^. 

The  following  is  a  general  property  of  reduction,  but  it  may  also  be 
checked  explicitly  in  this  case. 

3.3  Lemma.  If  F  and  G  are  functions  of  (m,v)  which  are  rotational! y 
invariant,  then 

{(F,S)>  =  =  const  C3.21! 

i.e.  the  full  KKS  bracket  is  the  same  as  if  it  were  computed  holding  m3 
constant. 

With  ilvli  =  1 ,  nvv  =  p,,  any  rotationally  invariant  function  of 
(m,v)  can  in  fact  be  expressed  as  a  function  of  V2IR1  -  v-jR^  and  v3> 

These  variables  are  a  convenient  representation  of  the  coadjoint  space 
reduced  by  the  S1  action.  Brackets  of  functions  of  these  variables  may  be 
computed  by  3.3.  Again,  one  can  check  by  hand  that  for  such  functions 

{{F,G}}  =  {F,G}0  (3.22) 

where 

(Fr\  -  3F  3G  9G  3F 

1  ’  '9  "  ae  ap  "  as  ap' 

y  9 
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§4.  Transversal  intersections  of  Invariant  manifolds  for  KKS 
variables  coupled  to  action  angle  variables 

We  now  develop  an  extension  of  the  Melnikov-Arnold  theory  (Kelnlkov 
[1963],  Arnold  [1964],  Holmes  [198(3)  which  applies  to  systems  described  by  a 
set  of  m  KKS  variables  u  =  (y-j ,  ....  u  )  £ and  a  set  of  n  action 
angle  variables  (0plj),  ...»  (6n»In).  This  theory  will  then  be  applied  to 
the  free  rigid  body  with  attachments  in  the  next  section.  For  the  heavy  top, 
a  more  sophisticated  version  is  needed  for  systems  whose  variables  do  not 
decompose  so  cleanly.  This  is  the  subject  of  Section  6. 

The  action  angle  variables  can  be  those  associated  with  motion  near  an 
elliptic  fixed  point  in  a  one  degree  of  freedom  system.  For  example,  our 
methods  apply  to  the  system  obtained  by  coupling  two  rigid  bodies  if  we  ex¬ 
amine  the  motion  near  a  homoclinic  orhit  in  one  (such  as  an  orbit  connecting 
(Q,2.,0)  to  (0,-2., 0)  in  Figure  1)  and  a  small  periodic  orbit  in  the  other 
(such  as  an  orbit  near  U,  0,0)  in  Figure  1).  The  coupled  system  then  will 
have  horseshoes  under  the  conditions  of  Theorem  4.3  below.  (We  believe  that 
this  applies,  in  particular,  to  the  five  mode  truncation  of  the  Euler  equa¬ 
tions  for  an  ideal  fluid  on  a  two  torus,  which  consists  of  two  sets  of  over¬ 
lapping  and  coupled  rigid  body  equations;  note  that  for  ty  »  so(3)  and  n  =  1 , 
we  have  five  variables  in  the  set  (p,e,I).)  For  simplicity,  however,  we 
have  chosen  a  rigid  body  with  attachments  to  work  out  in  detail  since  the 
action  angle  variables  are  more  explicit  and  it  is  easy  to  add  on  additional 
ones. 

We  assume  that  our  Himiltonian  takes  the  form 


H£(u,8,  I)  =  F(y)  +  l  G»(lj  +  eH*(u,e,  I) 

j»l  J  0 


=  Hu(y,  I 


In)  +  eM1 (p,6 ,1)  ♦  0(e2) 


(4.1) 


^  *  *  •  • » 


4.2 


* 

where  u  =  Gij.  .  ..t  Pm)  »  the  dual  space  of  a  Lie  algebra  and 

6  =  Cer  ....  en),  I  =  (1^,  •••*  In)  With  a  2Tr-periodic  variable. 

We  assume  that  the  KKS  system  associated  with  F  has  a  homoc! Inlc  (or 

_  * 

heterocllnlc)  orbit  y(t)  €<^  .  The  oscillator  frequencies 


3G . 


(4.2) 


are  assumed  to  be  positive. 


As  in  Holmes  and  Marsden  [1981btc]  we  can  solve  the  equation 


h  tu.e.i)  =  h 


for  I  In  the  form 
n 


1^  =  L  •••  >h)  *■'  L  Cu»  •  •  •  •  f 

+  eL  (y,  0^i  ...»  0^»  I i »  •••»  I^_^»h)  0(e  ) 

where 


Lu(y,Ir  ....  In.r  h)  =  6'J(h  -  F(y)  -  \  G^)) 

j  1 


L  (y»0-|»  •  ••»  *•** 

H  (y,0^j  •••»  0^,  •••»  t  (u»  ^i*  •••»  ^_i»h) 

nnU-0(y.  . In.r  h)) 


(4.3) 


(4.4) 


(4.5) 


(4.6) 


In  addition,  we  eliminate  t  in  favor  of  the  new  'time'  8n  and  write  ' 
for  d/den. 


4.3 


4.1  Proposition.  Hamilton's  equation  for  (4.1),  namely 
U1  =  {(yi,  H£}}  ,1-1,  .... 


A  _  3He  ;  .  3Hf  .  ,  _  , 

0j  "  3Ij  *  !j  36j  *  ^  *  ****  1 


(where  {{  ,  }}  denotes  the  KKS  bracket  In  the  u-variahles)  become 


(4.7) 


Uj  -  -«P,.  LE)) 

A '  =  -  — —  T  1  9  dl  ~ 

J  aI0  '  i  36J 


(4.8) 


Proof.  Implicit  differentiation  of  (4.3)  gives 


■SH£  +  3He  6L^_ 
5y 

It  follows  that  for  a  function  K.(y), 


=  0 


{{Me}}(y)  = 


,  r6K  6H%  . 

-< u-  is0  > 


.  r5K  3He  6Le -i  . 

s  [6?  n;^-1 > 


and  so  {{K,He)}/|r  =  -{{K,lG}}.  However  then,  K'  -  X  -  — ^ - -  -{{K,Le}> 
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6n  3HC/3I 


The  second  set  of  equations  in  (4.8)  follows  in  the  same  way. 


Using  l£  c  ♦  el^  +  0(e2}  we  see  that  (4.8)  has  the  form  of  a 
periodically  perturbed  system.  Next  we  relate  the  reduced  and  non-reduced 


brackets. 


4.2  Proposition.  We  have 


M-Vn  =  — {{f,h’}} 


14.9) 


Proof.  From  (4.5), 

3L°  „  1  6F 

6y  “  <5U 

and 


«n 


'  l  1  ol 
6H  f  3H  6LU 

.  3^  $U  , 


H 


1  3nn  6L^ 
6u 


However,  since  [ 


0 


or  6L_ 
5u  *  6u 


•J  =  0,  we  get 


,,0.1..  ,  ,«L°  SL1,  . 

(U  ,L  ))  -  ■< u.  L-^>  -jj-J  > 


a 


-<p. 


1 


nS  F  6H^  -i  . 
l6y’  6u  J  ' 


» 4  ■ 

Qn 


Let  us  now  give  a  special  case  of  the  general  result,  suitable  for  two 
degree  of  freedom  systems. 


4.3  Theorem.  Suppose  p(t)  is  a  hornoclinlc  (or  heteroclinic)  orbit  for 
F,  which  lies  on  a  coadjoint  orbit  in  V  of  dimension  2.  Furthermore. 
suppose  n  »  1 . 

Let  h  =  F(p)  be  the  energy  of  the  homoclinic  orbit  and  let  h  >  FT 
ajKl  i°  -  G_1(h  -  h)  be  constants.  Let  {{F.H1}}  (t,e°)  denote  the  KKS 
bracket  of  F(u)  and  Hl(yt  a(t°)t  ♦  5°,  t°)  evaluated  at  y(t).  Let 


4.5 


00 

Hie0)  - - V  {{F.HbKt.e0)  dt  (4.10) 

flnU0) 

and  assume  M(e^)  has  simple  zeros.  Then  for  e  >0  sufficiently  small. 
the  Hamiltonian  system  (4.1)  contains  transverse  homocllnlc  orbits  and  hence 
Smale  horseshoes  on  the  energy  surface  He  3  h. 


Proof.  By  reduction,  it  suffices  to  check  that  the  Melnikov  function  for  the 
reduced,  forced  system  on  the  coadjoint  orbit  containing  p  has  simple  zeros. 
This  involves  only  a  generalization  of  the  one  degree  of  freedom  Melnikov 
theory  for  forced  oscillations  to  two  dimensional  symplectic  manifolds.  The 
standard  proof  (see  Holmes  [1980]  and  Greenspan  and  Holmes  [1981]  )  carries 
over  directly.  The  Melnikov  function  for  the  reduced  system  is 


f  {{L°,l1»  de 


(4.11) 


since  ^  KKS  bracket  coincides  with  the  Poisson  bracket  on  coadjoint  orbits. 
Using  (4. hi  and  ~  52, 


Finally,  note  that 
theorem  results.  ■ 


H(8°)  - 


UF.H1}}  dt 


fl  °  Q(i^)  is  constant  on  the  homoclinic  orbit,  so  the 


To  deal  with  the  situation  in  which  n  >_  2,  we  introduce  the  following 
conditions  on  the  Hamiltonian  (4.1). 

(HI }  F  contains  a  homo-  (or  hetro)  clinic  orbit  d  with  energy  h. 

The  coadjoint  orbit  containing,  u  is  assumed  to  he  two-dirans tonal . 

The  saddle  points  for  y  are  denoted  (they  t-ould  be  coincident). 


4.6 


(H2)  QjUj)  =  Sj(Ij)  >  0,  j  -  \  ....  n. 

To  explain  the  remaining  conditions  some  discussion  is  needed. 


For  e  =  0,  note  that  the  Hamiltonian  system  for  has  two  n-^1 
parameter  families  of  invariant  n-1  dimensional  tori  T±(h1,  ...»  h^) 
given  fay 

P  =  P±  ) 

Gj(Ij)  =  hj  =  constant  (i.e.  Ij  =  =  6“^ (h^ ) )  \  (4.12) 

9j  =  flj^V6n  +  6j^  ^ 1180(1  2ir^’  ^  =  •••»  n-1  / 


(Correspondingly,  the  system  for  has  two  n  parameter  family  of  In¬ 
variant  tori  (T>(h-j.  . ...  hn)-)  Henceforth  we  write  the  (phase)  constants 
of  integration  8.(0)  as  0?,  j  »  1,  n-1,  n. 

J  J 

The  tori  T+(h1f  ....  hn  ^)  are  connected  by  the  n-dimensional  homo¬ 
clinic  manifold  defined  by 


U  •  u(9n  -  0n) 


GJ(Ij}  =  hj* 


(4.13) 


,0 


V  r‘j(Ij)en  +  9rj3l‘* 


...  n-1  , 


where  the  phase  constant  8^  associated  with  the  'reduced*  degree  of  freedom 

appears  explicitly.  This  manifold  consists  of  the  coincident  stable  and 

unstable  manifolds  of  the  tori  T  (h. ,  ...»  h  ,);  i.e. 

i  1  n- 1 


WS(Tt(h1 . hnl))  *  MU(T,(h1 . hnl)) 


given  by  (4.13). 


4.7 


For  e  ^  0  the  system  (4.7)  possesses  a  Poincare  map  Pg  from  (a 
piece  of)  (p,  0^,  ...,  0^,  1^,  ....  I  -j)  space  to  Itself  where  0n 
goes  through  an  Increment  of  2it,  starting  at  some  fixed  value  0®,  (which 
will  be  suppressed  In  the  notation).  Below,  when  we  refer  to  transverse 
intersection  of  stable  and  unstable  manifolds,  we  mean  so  for  this  Poincare 
map. 

(H3)  Assume  that  the  constants  G . ( I . )  *  h.,  j  =  1,  ....  n  are  chosen 
so  that  the  unperturbed  frequencies  ilj(I^),  ....  &|(In)  satisfy 
the  non-degeneracy  conditions  (i.e.  Ql(I.)  f  0,  j  =  1,  ....  n-1) 

J  J 

and  the  non-resonance  conditions  of  the  KAM  theorem,  (cf.  Arnold 
(1978),  Appendix  8). 


This  condition  ensures  that  the  tori  T+(h^,  ...  h^  ^)  perturb  to 
invariant  tori  T  +(h.,  ....  h  ,)  for  for  e  sufficiently  small. 

_  n 

Let  h  >  R-,  h  3  h  +  £  h.  where  h.  >  Q,  and  the  unperturbed  homoclinic  mani> 

j=l  J  J 

fold  be  filled  with  an  n-parameter  family  of  orbits  given  by  (u,  9^,  _ _ 

9n,  Ir  In)  =  u(t) ,  Qj ( Ij  )t  ♦  9®,  ...,  «nUn)t  +  0°.  Ir  ....  In). 

Pick  one  such  orbit  and  let  denote  the  KXS  Poisson  bracket  of 

F(y)  and  H^(y»  8^,  ....  0^,  1^,  ....  1^)  evaluated  on  this  orbit.  Sim¬ 
ilarly,  let  U  .H1}  =  ,  k  5  1,  ....  n-1  be  evaluated  on  this  orbit. 

Define  the  Heinikov  Vector  N(88)  .  (H^ ,  .....  M  H^)  by 


V8? 


Sn‘  h’  V  h2 . Vl}  a  j  {Ik*  *]}  k  *  1.  n-1 


Mn(6l*  *•••  9n*  h‘  V  h2 . h  )  »  4“  J* 


HF.H1}}  dt 


(4.14) 


4.8 


(We  note  that  h. 


n-1 

h  -  -  l  h 


,,  I 
J  n 


and  hn  do  not  explicitly  enter 


the  calcualtlons,  since  In  Is  eliminated  by  the  reduction  process;  we  also 
note  that  these  Integrals  need  not  be  absolutely  convergent,  but  we  do  require 
conditional  convergence.) 


(H4)  Assume  that  the  multiply  2tt  periodic  Melnikov  vector  M  :  IRn  ♦  IRn 
(which  Is  Independent  of  e)  has  at  least  one  transversal  zero;  I.e.  there 
is  a  point  (0°,  ...,  0°)  for  which 

M(9?,  ....  0°)  =  0 

—  det[DM(0°,  ....  0°)]  f  0  ,  ' 

where  DM  is  the  n  *  n  matrix  of  partial  derivatives  of  Hj,  ....  H  with 
respect  to  0°,  .....  0",  the  initial  phases  of  the  orbit. 


Here  is  the  result  for  n  >  2. 

4.4  Theorem.  If  conditions  (HI)  -  (H4)  hold  for  the  system  (4.1),  then, 
for  e  sufficiently  small,  the  perturbed  stable  and  unstable  manifolds 

WS(T  ),  and  WU(T  of  the  perturbed  tori  T  intersect  transversely. 

e.t  •  1  1  £ r  "  " ■  ■ r ■ . '  1  "r ■  c ,t  1  1 

The  theorem  follows  from  the  arguments  of  Holmes  and  Karsden  [1981c] 
in  the  present  context.  We  also  refer  the  reader  to  that  paoer  for  a  discus¬ 
sion  of  how  this  yields  Arnold  diffusion  and  for  related  references. 

There  is  a  similar  result  when  the  coadjoint  orbit  is  higher  dimensional 

* 

i.e.  of  dimension  2H,  H  >  1,  but  the  system  for  F  on  a  is  completely 


4.9 


Integrable,  say  with  integrals 

F  a  Fj  ,  f2t  Fm  « 

where  ....  FM  are  associated  with  action  angle  variables.  Now  a 
result  similar  to  4.4  holds  If  the  Melnikov  vector  Is  enlarged  by  replacing 
the  number 

jf  {{F,H*}}  dt 

“n  '~ca 

by  the  vector 

00 

J-  {{F^H1}}  dt,  1  -  1,  ....  m 

n 

ti' would  be  of  interest  to  apply  such  a  generalization  to  the  Toda  lattice 
and  related  completely  integrate  systems. 


IS.  An  example:  the  rigid  body  with  attachments 


We  start  by  considering  the  rigid  body  in  the  absence  of  gravity,  with 
a  single  attachment  which  spins  freely  about  an  axis  coincident  with  one  of 
the  body's  principal  axes,  say  1.  The  angle  about  this  axis  is  denoted  e. 
For  simplicity  we  suppose  that  the  attachement  does  not  affect  the  position 
of  the  principal  axes,  so  that  the  inertia  tensor  remains  diagonal,  but 
that  the  attachment  is  slightly  asynmetrical  in  the  sense  that  the  perturbed 
moments  of  inertia  about  the  2  and  3  axes  depend  upon  the  angular  position 
8  of  the  attachment  relative  to  the  body.  Such  an  attachment  could  consist 
of  a  (heavy)  disc  and  a  light  rod,  as  in  Figure  4. 


2 


Figure  4 


Letting  the  moments  of  inertia  of  the  body  with  respect  to  the  three  body 
axes  be  Ij  j  =  1,2,3,  and  of  the  attachment  be  0^  +  e.  Or,  +  c  cos^e, 

♦  e  sin  8.  We  obtain  the  Hamiltonian 


5.2 


ue  -  1 
H  -  o 


1 

\  -5 _ + 

7  i!  +  e 


"2  m2  I2 

I.  +  c  cos29  I,  +  e  s1nz«  J1  +  e 


7  j  «j/ij *  i2/2J,  -  f  fe  A  ^  s1"2e  *4j *  ou2) 

Ja|  l1!  x2  *3  dl  J 


=  F  +  G  +  eH  +,0(.e2) 


(5.1) 


where  I.  =  U  +  <3.  and  where  (8,1)  are  action  angle  variables  for  the 
J  J  J 

2  1 

attachment.  Since  the  unperturbed  system  is  a  produce  flow  on  (S  )  x  (ir  x  s  ), 

the  product  of  the  coadjoint  orbit  of  S0(3)  with  the  (1,8)  cylinder,  we  can 

use  (3.9)  to  write  the  homoclinic  orbits  for  an  energy  level  =  h  a 
(2  ?1 

1  cf  ir-  2  2  2  2 

2  T~  +  J-  *  where  ^  3  ^  +  +  m3  and  k  is  a  constant,  as 

l12 

/Sj 

m-j  =  ±i  ~~  sech(-va^a3il  t) 

^  2 


mg  =  ±t  tanh(-^ara3t  t) 

"3  =  ±lM^  t) 


(5.2) 


I  =  k(constant) 


9  =  ~  t  +  8° 
J1 


To  show  that  transverse  homoclinic  orbits  occur  for  &  f  0  we  need 
only  show  that  the  Melnikov  function 


M(90>  ■  W) 


{{F,H‘}}  dt 


(5.3) 


has  simple  zeros,  all  other  conditions  of  Theorem  4.3  being  immediately 


satisfied.  Note  that  we  must  set  k  >  0,  so  that  I  >  0  for  the  unper¬ 
turbed  system  and  hence  the  inversion  of  He  *  h  goes  through.  This  is 
the  reason  for  choosing  an  attachment  with  a  heavy  'in  balance*  symmetric 
component).  The  KKS  bracket  is  given  by  (3.6): 


{{F(m),  HVl>  6-6°) }}  -  -m.ymF  x  ^H1 


a^  a2  cos2©  a^  sinks'1 

~l  +  — 7 - +  — 


(5.4) 


Dlj 


Noting  that  on  any  homoclinic  orhit  and  m3  are  even  while  m2  is 

2 

odd,  it  follows  that  the  constant  term  (aj/I.|)m.m2m3  vanishes  in  the  Melnikov 
integral  and  we  are  left  with 


M(e°) 


1 


f 


W7T  Umlm2m3 

'a 

mlm2m3  ~ 

2 


~i  cos20  +  sin2e 
*2  *3 


dO, 


2k  jmW  mlm2m3 


a«  a« 

—k  cos  2s  -  — k  cos  20 

1°  I3 


do. 


_L 

2k 


f*2 

T  ‘  T 

1*2  *3J 

. 

m^m2m3  cos  20  d8. 


(5.5) 


Inserting  the  expressions  for  the  homoclinic  orbits,  we  get 


M(e°)  .  ^ 


/ 

00 

<Zr— -) 

h. 

( 

li 

a3j 

i  -00 

i  -*2  j 

sechz(-^a"3's,  t)  tanh(-v^a^t  t) 


cos  2 


37t+e° 


dt  , 


?  _  _  fOL.f'S 

=  C!  sech^-v^a^  t)  tanh(-v^ajit  t)  sinj^pj  dt 


sin  20°  ,  (5.6) 


5.4 


where 


°1 

2iT 


a2  a3 
- 7 

U2  *1 


The  integral  of  (5.6)  may  be  evaluated  by  the  method  of  residues  to  obtain 


M(0°)  =  C  —-X — - t  cosech 

j^(-y£^ajir 


kir 


-*V35'  J1 


sin  26°  ,  (5.7) 


which  has  simple  zeros.  We  therefore  have 


5.1  Theorem.  The  free  rigid  body  with  a  single  slightly  asymmetrical 
freely  rotating  attachment  on  one  of  its  principal  axes  possesses  transverse 
heteroclinic  orbits  and  hence  Smale  horseshoes  in  a  suitably  chosen  cross 
section  of  the  constant  energy  surface  with  k  >  0. 


This  implies  that  the  rigid  body  equations  with  an  additional  attach¬ 
ment  are  non  integrate.  More  precisely,  if  we  make  a  Markov  partition  of 
the  invariant  sphere  consisting  of  the  four  open  regions  filled  with  periodic 
notions  in  the  unperturbed  case  (Figure  1),  then  the  dynamics  of  the  per¬ 
turbed  Poincare  map  is  conjugate  to  the  subshift  of  finite  type  on  these  four 
symbols.  To  see  this  we  sketch  the  homoclinic  structure  on  the  sphere  in 
Figure  5,  identifying  one  of  the  centers  (in  region  A)  with  the  point  at 
infinity.  It  is  clear  that  orbits  starting  near  the  manifold  on  the  'boundary* 
of  regions  2  and  3  can  be  selected  such  that  they  pass  either  from  region 
2  +  2  or  2  +  3  or  3  +  2  or  3  +  3.  Similarly  on  the  border  of  1,2  orbits 
can  be  found  passing  from  1  +  1,  1  +  2,  2  +  1  or  2  +  2.  Continuing  in  this 


Figure  5.  The  homoclinic  structure  of  the  perturhed  manifold 
^(O.ii.Q),  Wn(.0,±2.,0)  on  the  sphere. 

way  we  find  the  transition  matrix 

1  1  0  1  ' 

1110 
0  111 
10  11 

where  a^  -  1  if  there  is  an  orbit  from  region  1  to  region  j  and  a^-  a  0 
if  there  is  no  such  orbit. 

If  two  or  more  attachments  are  added  to  the  free  rigid  body  then 
Arnold  diffusion  may  take  place.  Taking  a  system  with  two  identical  attach¬ 
ments  on  axis  1,  each  free  to  rotate  independently,  with  momenta  1,0  and 
angles  9,  <p,  we  have  the  Hamiltonian 


A  =  = 


5.6 


=  F  +  G  +  G  +  eH1  +  0(e2)  -iff 

'  1  U=1  J1 

■  f  l"7+  “§  Ccos2e  +  C0S24>)  +  Csin20  +  s1n2«J>)  +  - — ]  +  0(e2) 

I1?  *2 


(5.8) 


In  this  case  the  Melnikov  vector  consists  of  the  pair  of  function 


GO 

M1(e°,«t>°)  =  [  (i.Hbdt 

'  tmCO 

00 

¥6°-*°>  =Lwum’,k 


1  3H 1  J 

where  {I,H  }=  -  -gg-  is  the  usual  canonical  Poisson  bracket  and  «2  =  . 

Computations  similar  to  those  above  (also  cf.  Holmes  and  Marsden  [1981c],  §4) 

show  that 


M-j  =  C-|  sin  2e1 

M2  =  C2  sin  20°  +  Oj  sin  20° 


(5.10) 


where  the  C,  are  non-zero  constants  depending  on  I,  J,  Ji,  and  the  choice 

J 

2  2  2 

of  total  energy  and  energy  in  each  'mode'  (E  ■  S.  /2I2  +  k^/2J^  +  k2/2J^). 

Thus  M  has  simple  (transverse)  zeros  for  0°  »  nm,  <J>°  «  mr,  m,n  G  2Z  . 

Since  Ql(I)  =  —■  t  0,  all  the  conditions  (HI)  -  (H4)  of  Theorem  4.4  hold. 

1  J1 

Therefore,  the  rigid  body  with  two  attachments  as  specified  exhibits  Arnold 
diffusion.  The  existence  of  a  transition  chain  of  two- tori  connected  by 


5.7 


heteroclinic  orbits,  (see  Holmes  and  Marsden  [1981c],  S3  for  a  discussion) 
Implies  that  angular  momentum  can  be  transferred  back  and  forth  between 
the  two  spinning  attachments  in  a  chaotic  manner. 

Remark.  An  amusing  corollary  for  the  case  of  a  single  attachment  Is  that 
there  are  configurations  of  the  Euler  elastica  for  which  the  sequence  of 
loops  above  and  helow  the  mean  level  can  be  prescribed  in  advance.  (For 
example  the  loops  can  be  coded  by  the  binary  expansion  of  an  irrational 
number).  This  follows  from  the  above  calculations  and  the  remarkable  fact 
that  the  elastica  equations  have  the  form  of  the  equations  of  a  rigid  body 
with  an  attachment;  see  Love  [1927,  p.  400]. 


6.1 


§6.  Melnikov  Theory  for  Systems  with  Symmetry 

* 

Me  now  develop  a  version  of  the  Melnikov  theory  that  applies  to  per- 
turbations  of  a  two  degree  of  freedom  system  with  an  S  symmetry.  We  have 
chosen  this  context  with  applications  to  the  motion  of  a  nearly  symmetric 
heavy  top  in  mind. 

The  key  new  feature  is  that  the  unperturbed  system  is  no  longer 
assumed  to  be  a  product  system  consisting  of  variables  with  a  homoclinic 
orbit  and  action  angle  variables.  Rather,  this  product  structure  is  gen¬ 
eralized  to  the  assumption  of  an  S1  reduction. 

Roughly  speaking,  our  umperturbed  Hamiltonian  HQ  no  longer  can  be 
split  as 

H°(q,p,D  =  FCq.p)  +  G(.I) 

so  that  the  frequency  function  Q  s  now  may  depend  on  (p,q).  This 
is,  in  fact  the  situation  for  the  nearly  symmetric  heavy  top. 

Let  us  start  with  a  four  dimensional  symplectic  manifold  P,  whose 
points  are  denoted  x.  Suppose  acts  on  P  by  canonical  transformations 

if 

and  has  an  Ad  -equi variant  momentum  map  J  :P  IR.  Let  the  reduced  space 
be  denoted 

Pu  = 

(see  Marsden  and  Weinstein  [1974]). 

*  2  1 

For  the  heavy  top,  P  *  T  S  and  S  consists  of  rotations  about  the 

axis  of  symmetry.  To  keep  the  notation  consistent,  we  shall  use  y  for  the 
1  ] 

angle  on  S  .  Motivated  by  this  example,  we  allow  Py  to  have  isolated 
singularities,  but  in  this  case  we  demand  that  the  constructions  carried  out 
i  below  make  sense  at  the  singular  points.  For  the  heavy  top  this  causes  no 

f: 

r  difficulties. 

E 

f 

& 

?■ 


6.2 


Let  points  in  the  reduced  space  be  denoted  u  E  P  .  Thus,  u  consists 
of  an  S  orbit  in  P,  lying  in  the  level  set  0-1(y).  Choosing  a  slice 
(cross-section)  for  this  action,  parametrizes  the  point  on  the  orbit  and 
of  course  the  value  of  0  labels  the  surface  J"\y).  Thus,  we  write  points 
x  as 

x  =  (u.tM)  (6.1) 

See  Figure  6. 

In  this  notation,  a  function  of  x  is  invariant  if  and  only  if 

it  can  be  expressed  as  a  function  of  u  and  J  alone.  Such  functions  induce 

functions  on  P  ,  the  reduced  space,  and  have  well  defined  Poisson  brackets 

on  p  .  In  fact,  the  Poisson  bracket  of  two  such  functions  is  just  that 

function  induced  on  P  hy  their  Poisson  bracket  on  P.  Of  course  any  S^ 

V 

invariant  Hamitonian  on  P  gives  a  completely  integrable  system,  the  integrals 
being  J  and  H,  or  equivalently  H  for  the  reduced  one  degree  of  freedom 
system. 


1 


Figure  6 


6.3 


How  assume  that  we  start  with  an  Invariant  Hamiltonian  H^(u*J) 
and  we  perturb  it  by  a  general  Hamiltonian.  Thus,  write 

He(u,iM)  =  H°(u,J)  +  eH1Cu,«,J)  +  0(e2)  C6.2) 

Now  we  must  generalize  the  procedure  of  Holmes  and  Marsden  [1981b]  to  this 
context. 

We  begin  by  letting 

flCii.O)  =  -gj  (u,J)  (6.3) 

and  assuming  a  is  positive,  so  that  H°(u,0)  is  invertible  In  the  J- 
variable.  For  e  small  then 

He(u,^,J)  »  h  [6.4) 

can  he  solved  for  J.  Write  H^|(J)  a  H°(u,J)  and  (H^)'1  for  Its  inverse 
in  the  J  variable. 


6*^  Proposition.  J  =>  Le(u,^,h)  =  Lu(u,h)  +  el  (u,<j>,h)  +  0(e£)  where 


and 


L°(u.h)  ■  (Hj)*’(h) 


n(u»l  (u,h)) 


(6.5) 


(6.6) 


Proof.  (6.4)  reads 


He(u,*,  1°  +  eL1  +  0(e2})  =  h 


i  .e. 


H°(u,  1°  +  dJ  ♦  0(e2))  +  zh\ u,  i|»,  l^  +  eJ  +  0(e2))  =  h 


6.4 


i.e.  H^(u,L°)  +  ca(u,t°)  +  eH^u.^.L0)  =  h  +  0(e2).  Comparing  powers  of 
e  gives  the  result.  M 


As  in  Marsden  [1981,  Lecture  4]  we  can  arrange  things  so  that  tp  and 
J  are  conjugate  variables.  Thus  under  the  dynamics  of  He, 


.  MT 


(6.7) 


♦  =  4.H  }  -  -jjj- 

Now  we  change  time  variables  from  t  to  f  and  write  ’  for  d/dip.  Thus, 
for  a  function  F(u), 


F‘  =  ~ 

However,  from  (6.4)  and  properties  of  the  Poisson  bracket, 

(F(u),  He(u,ip,L~(u,ip»h) ) >  =  0 


so 


{F,He}  +  {F,Le}  =  0 


Comparing  (6.8)  and  (6.9), 


(6.8) 


(6.9) 


F*  =  -  ( F,Le }  »  -{F,l°}  -  eiF.L1)  +  0(e2) 


(6.10) 


Thus,  (6.10)  is  in  the  form  of  a  periodical 1>  forced  Hamiltonian  system  on 
the  reduced  phase  space.  The  brackets  in  (6.10)  are  taken  in  the  u-variable 
alone.  Thus,  we  have: 


6.2  Proposition.  The  evolution  of  u  as  a  function  of  is  a  ^-dependent 
Hamiltonian  system  in  the  reduced  phase  space .  The  equations  of  evolution 
are  given  by  (6.10). 

Now  suppose  that  the  reduced  system  for  e  =  0  h&s  a  homoclinic  orbit. 
To  detect  the  presence  of  transverse  nomoclinic  orbits  for  e  0  we  must 


* 


integrate  around  this  orbit.  Now  in  examples,  and  0  and 

their  bracket  can  be  laborious  to  compute.  It  will  save  some  effort  if  we 
relate  to  H°  and  H*. 

6.3  Proposition.  (L0,^ }  =  ^  (H0,  (6.11) 

where  ft  is  evaluated  at  u,  J  and  the  bracket:  are  taken  with  respect  to 
u,  holding  (j/,  J  fixed. 

Remarks  1.  Formula  (6.11)  is  a  generalization  of  (3.1)  of  Holmes  and  Marsden 
[1981b].  In  the  later  case  ft  was  independent  of  u. 

2.  The  above  development  does  not  depend  on  the  dimension  of  P. 


Proof  of  6.3.  Given  F(u),  we  first  compute  (l^,F}.  From  (6.5), 
H°(u,L°(u,h))  -  h,  so 

{H°,F}  +  ftU.°,F}  -  0 


Thus 


{L°,F}  =  -  i  {H°.F> 


(6.12) 


Next,  from  (6.6)  we  obtain 


(r  ,1,  H1,  1  3H1  {r  .0,  .  H'  3ft  ,0. 

{G,L  }  -  -{G,  — }  -  q  -jj  (6.1  )  +  jr  (S.L  > 


,1 


(6  13) 


Taking  G  s  in  (6.13]  the  last  two  terms  vanish,  so 

,.0.1,.  0  H], 

(L  ,L  )  s  -{L  , 


By  (6.12)  this  gives  the  stated  result. 


Assentling  these  results  as  we  did  in  Section  4,  we  obtain  the  following. 


6.6 


6.4  Theorem.  Consider  a  two  degree  of  freedom  Hamiltonian  system  on  a 
symplectic  manifold  P  with  an  symmetry  with  a  Hamiltonian  of  the 
form  (6.2).  Assume  that  J  is  chosen  and  fixed  so  that  the  reduced  system 
H^(u»J)  has  a  homoclinic  (or  heteroclinic)  orbit  u(t)  in  the  reduced  space 
Pu  and  so  0(u(t)#J)  >  Q.  Let 


4»(t)  3 


ft 

•0 


Q(u(t),J)  dt  + 


(6.14) 


0  H 1  0  _ 

and  let  {H  ,  — }  (t,4>  )  denote  the  u-Polsson  bracket  evaluated  at  u(t), 

<t»(t)  and  J.  Let 


(H°,  i~)  (t,*0)  dt 


(6.15) 


and  assume  M(^)  has  simple  zeros  as  a  function  of  Then  for  e 

sufficiently  small,  the  system  (6.2)  has  transverse  homoclinic  (or  hetero¬ 
clinic)  orbits,  and  hence  Smale  horseshoes  on  the  energy  surface  He  *  h, 
where  h  c  H°(u,J), 


7.1 


.'•v 


§7.  Example:  the  motion  of  a  nearly  symmetric  heavy  top 
We  now  show  that  Theorem  6.4  Implies: 

7.1  Theorem.  If  I^/I^  is  sufficiently  large,  1^  *  I-j  +  £  and  e  is 
sufficiently  small,  e  t  0,  then  the  Hamiltonian  system  for  heavy  too 
(see  3.11  and  3.12)  has  transverse  homoclinic  orbits  in  the  Poincard  map 
for  the  <p  variable  on  energy  surfaces  close  to  the  homoclinic  orbit 
described  in  3.2  or  3.4. 


7.2  Corollary.  The  heavy  top  close  to  the  symnetric  top  has  no  analytic 
integrals  other  than  the  energy  and  angular  momentum  about  the  vertical 
axis. 


Remarks  1.  As  we  have  already  discussed,  this  corollary  has  recently  been 

obtained  oy  Ziglin,  but  by  rather  different  methods.  Moreover,  our  result 

7.1  shows  the  existence  of  'chaotic1  orbits. 

2.  I  /I-  being  large  can  be  replaced  by  0  <y  <1  and  the 
1 

integral  7.13  below  being  non-zero.  This  integral  is  non-zero  for  most 
values  of  M,  6,  1^,  Iv  as  we  shall  show. 


lie  shall  prove  theorem  7.1  In  the  Euler  angle  representation  first 
and  then  sketch  how  the  proof  can  be  alternatively  obtained  using  the  KKS 
description. 

In  (3.12)  let  I^  a  1^  +  e  .  This  gives 

,  2  ?  2 
0  1  UP* '  P*  cos  ^  Pq  p  J 

H  £  I  ]—_■  T- -  *  -  +  T-!  +  Hg4  cos  0  (7.1) 


*1  sin  9 


‘1 


■IS 


7.2 


and 

H1  - J — y-  ((p.  -  prf  cos  6)cos  if  -  p  sin  0  sin  ^)2  (7.2) 

2I|  sl/Te  *  ^  6 

Note  that  Is  the  Hamiltonian  for  the  symmetric  top  and  so  has  a  homo¬ 


clinic  orbit  given  by  (3.15). 


Since  <j»  and  p.  are  conjugate  variables  and  J  =  p. ,  we  have  from 

y  V 


(6.3), 


a  « ii-  - 

a3 


Ij  siA 


*  cos  0 


(7.3) 


On  the  homoclinic  orbit,  »  constant  =  I^b,  so  (7.3)  becomes 


a  =  b 


1  cos  0 


I3  1  +  COS  0 


(7.4) 


3  Lemma.  With  p.  *  p,  =  I,b,  8  *  ^-v^—  , 

- -  4>  ip  i  ii 


♦  *  T'+^osT  *21  6  *  if  * 


we  have 


(H°*  ir}  s  ~  (A(e)  sin  2^  +  B(0)  cos  2<p  +  C(e) )  (7.5) 


9»Pc 


where 


2:2, 


b  (12  8  . 


A(e)  *  f  ir  -  * + 


2  ^  ’  2  '  bn 


(1  -  COS  0)  -  0 


.2 


B  ( Q ) 


4  *  ®  *<«  •  - & 


and 


C(e)  a  0 


•4  .  >2 -2) 

§. .  »s.-y*ec<r 

4  4b£ 


sin  0 


•  \  ’ 


.1 


Proof,  The  computations  are  slightly  tedious  hut  straightforward.  He 
write 


Ik°,  4U  mV}  -  —  in0 


Q 


or, ay  , 


rH°  Hl,  _  9l£  3h1  SH°  3H1 
1  *  i  ~  36  3pQ  *  3p  36  ’ 

/H°  a\  - 

{H  >n}  ~  “  3pQ  3e  9 


and  compute  that 


3H°  _  f*2  .  A  6  •  i 

-gjj*  *  -Ii  sm  0  »  ^  SIR  6  , 


MI 

3P, 


=  § 


(7.6) 


M_  =  ^  sl-n  9  cos2,{,  +  &  Sin  ?<!», 

=  ~  (I  sin  6  sin  Zip  -  0sin2\j>), 
3Pe  ^ 


|  [($2  -  |)  sin2e  -  02]  sin  2i|> 


and 


+  e^sin  e  cos  Zip  -  M  sl-n  e(cos  2ip  -  1), 

to  -  £  cin  fl 

W  "  T  sin  6  ’ 

•  -2 

{H0,^}  =  -  sin  6. 


Substituting  these  expressions  into  (7.6)  and  simplifying  yields  (7.5). 
The  Melnikov  function  is  given  by 


M(ip°) 


K  4! « 


0-7) 


7.4 


where  the  integral  is  evaluated  along  the  homoclinic  orbit 
cos  0  ®  1  -  y  sech^  ■  ,  pQ  *  I-j§ 

Csee  (3.15))  and  where 


(7.8) 


r  ^  - 

Kt)  =  n(t)  dt  +  r  a  ?(t)  +  r 
jo 


(7.9) 


and 


fl(t)  *  b 


1  -  ysech^ 

[m 

Y  sech^ 

m  1 

(7.10) 


Note  that  if  0  <  y  <  1.  or  if  Ij/I^  is  sufficiently  large,  Q(t)  >  0. 
Substitute  (7.5)  and  (7.0)  into  (7.7)  yields 


L'  -oo 


(  i-CA(Q)  sin  2ty+B(0)  cos  2$)  dt 

oo 

1 


cos  2ip 


.0 


^  (A(0)  cos  2\p  -  B( 0)  sin  2\p)  dt 


sin  2<i>  + 


1 


j  a 

'  mOO 


kCO)  dt 


(7.11) 


The  first,  second  and  fifth  terms  are  odd  functions  of  t  and  so  we  obtain 
the  following 


7.4  Lemma. 


MOP0} 


■[i 


jj  (A(e)  cos  2i|> 


-  B(0)  sin  2ip)  dt 


;in  2iJ> 


(7.12) 


7.5 


.a  U 

Now  we  observe  that  at  t  =  0,  A  a  >  0  and  S  ■  0.  It  follows 

that  for  I-j/I3  sufficiently  large,  the  portion  of  the  integral  from  the 
first  term  near  t  =  0  dominates  and  so  the  integral  must  be  a  non-zero 
number.  Thus  M(^)  has  simple  zeros  and  the  theorem  is  proved. 


h2  ? 

Remark.  For  any  given  b,  6,  Ip  I3  (with  y  c  2  -  ~  and  b  <  23, 
b,  8  >  0)  one  needs  only  the  condition  that 

00 

r  1  —  _ 

jj  (A(6)  cos  2}p  -  B(6)  sin  2i|>)  dt  (7.13) 

—OO 


be  non-zero.  Since  we  cannot  evaluate  (7.13)  analytically  it  does  not 
seem  so  simple  to  decide  exactly  when  (7.13)  vanishes.  Since  we  can  be  sure 
it  is  non-zero  for  Ij/I3  large,  it  follows  that  (7.13)  can  vanish  for  at 
most  a  finite  set  of  values  of  b,  3,  I-j ,  I3  as  it  is  analytic.  Thus,  we 
can  be  sure  of  transverse  hmoclinic  orbits  for  generic  M,  Ip  I3,  b,  if 
0  <  y  <  1 . 

Finally,  we  indicate  how  the  same  computations  can  be  done  using  the 
KKS  formalism.  This  actually  makes  the  computations  slightly  easier,  but  the 
final  result  is  the  same. 

,  3  m2 

Again,  letting  =  I,  +  e  in  H  =  7  T  J-  +  Mga,  v~,  we  have 
2  2  J=1  3 


2  2  2 

n  m1  +  m?  m, 

Hu  =  -l-r—  +  of-  +  Mg*  v- 

1  a3  J 


21 


(7.14) 


and 


m 


2 


(7.15) 


7.6 


Next,  observe  that  tp  and  m3  are  conjugate  variables,  where  tan  ^  *  v^/Vg, 


so  If  m*v  s  I^b  =  m^, 


a  =  {{^,H0}}  =  b 


fIl  v3 


*3  ^  +  v3 


,0  H' 


(7.16) 


which  agrees  with  (7.4).  To  compute  {{H  ,  “}}  we  write  (see  Lemma  3.3) 

H1  as  a  function  of  the  reduced  variables  Vgtn^  -  v^g  and  v3  and  compute 
the  bracket  holding  i|>  and  m^  fixed.  This  is  done  by  writing 


2  1  ? 

m2  =  (m«v  "  ra3V3^  cos  ^  ’  (-v2ral  ■v]ra2^m’v  "m3v3^  s^n  ^ 

1  "v3 

+  (vgm-|  -  v^)2  sin2^ 


(7.17) 


When  nrv  =  =  I^b,  this  becomes 


m 


\  «  I,b  .  i1b(l  -  VjKvjin,  -  v,^)  sin  ,  +  (v^  -v,m2)‘  s1n‘i|i 

CM8) 

Q 

Using  these  expressions  we  readily  compute  {{H  ,  — }}  using  (3.17).  After 
substituting  from  Table  1,  the  same  expression  (7.11)  as  above  results. 

Thus,  the  proof  may  be  completed  in  the  same  way  as  with  the  Euler  angle 
argument. 


i2  _ .  2 , 


« 
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